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The vector and pseudoscalar decay constants are calculated in the framework of the Field Correla- 
tor Method. Di-electron widths: T ee (J/ip) = 5.41 keV, r ee (i/>'(3686)) = 2.47 keV, r ee (V"(3770)) = 
0.248 keV, in good agreement with experiment, are obtained with the same coupling, a s — 0.165, in 
QCD radiative corrections. We show that the larger a s = 0.191 ±0.004 is needed to reach agreement 
with experiment for T 77 (»7 C ) = 7.22 keV, F 77 (x( 3 -Po)) = 3.3 keV, F 77 (x( 3 P2)) = 0.54 keV, and also 
for T(J/ip -» 3g) = 59.5 keV, Y(J/ip -* j2g) = 5.7 keV. Meanwhile even larger a s = 0.238 gives 
rise to good description of T(ij)' — > 3g) = 52.7 keV, T(ip' — > ^/2g) = 3.5 keV, and provides correct 

3(j/V>-dight hadrons) _ g 24 



ratio of the branching fractions: 



I. INTRODUCTION 



B(*'->light hadrons) 



Low-lying states of heavy quarkonia have been an 
important laboratory to study both perturbative and 
nonperturbative phenomena in QCD. However, recent 
discoveries of higher resonances, in particular A (3872) 
[lj and y(4260) [1J have shown that these and some 
other new resonances cannot be interpreted as conven- 
tional QQ mesons. To understand the nature of new res- 
onances, evidently, two- (or many-) channel consideration 
is needed. However, in strict sense it cannot be done now 
because nonperturbative theory of strong decays is not 
still well developed in QCD. Therefore for identification 
of new resonances with J PC = 1 , observed in e + e _ via 
the initial state radiation Q, a special role belongs 
to di-electron widths and also two-photon widths for C- 
even resonances, which are reasonably well described by 
existing QCD formulas. At this point it is worthwhile to 
remind that the di-electron width of a QQ meson is by 
two orders (may be even more) larger than di-electron 
width of a compact four-quark system [f|. 

In our paper, firstly, we calculate the decay constants of 
vector (V) and pseudoscalar (P) mesons in charmonium 
using the Field Correlator Method (FCM), which has 
been successfully applied to heavy- light mesons @. Due 
to relativistic corrections di-electon widths and their ra- 
tios, calculated here, agree with experiment with high ac- 
curacy. Therefore from the absolute values of di-electron 
widths some important factors, containing the squared 
wave functions at the origin, can be extracted and then 
used in different annihilation decays. 

We pay a special attention to the influence of radiative 
corrections on different annihilation rates. The absolute 
values of di-electron widths are shown to agree with ex- 
perimental numbers only if the QCD radiative corrections 
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are taken into account. Unfortunately, at present there is 
no consensus about the true value of the strong coupling 
in them. These corrections, known in first (one-loop) ap- 
proximation > 9 > en ter the di-electron and two-photon 
widths as separate factors: (3y = 1 — , [3p = 1 — 3J " 3 . 
In [9(] these factors are put equal unity: (3y = f3p = 1.0, 
i.e. the QCD correction is neglected, while in [Tol ]- 
[i~3| their values are almost two times smaller, (3y = 
0.52 ± 0.06. The reason of this uncertainty partly oc- 
curs because the contribution of higher corrections re- 
mains unknown. Therefore, although by derivation the 
coupling a s in different annihilation widths is defined at 
the standard scale fi = 2mg or fx = M v r P \ [l4| (in the 
MS scheme), factually, this strong coupling appears to 
be an effective one and can differ in different annihilation 
decays, since for them higher order perturbative correc- 
tions can be different. 

In our paper we show that in the ip- family the di- 
electron widths are described with the same coupling, 
which turns out to be relatively small: a s = 0.165 
or (3v = 0.72 (the same "universality" is observed in 
bottomonium [la]). Meanwhile to describe two-photon 
widths of rye, Xc(l 3 -Po)j Xc(l 3 ^2), and also three-gluon 
annihilation rate of J /ip only the choice of larger cou- 
pling, a s — 0.191 ± 0.004, gives rise to agreement with 
experiment. Even larger a s = 0.25(2) provides correct 
number for the ip 1 width T(ip' — > 3g). Thus our analy- 
sis shows that low-lying charmonium states have no an 
universal scale for different annihilation decays and there- 
fore any ratio of their widths cannot be used to extract 
characteristic strong coupling (for the discussion see [l6| , 
[13]); m particular, they are different for the J/ip and ip' 
three-gluon annihilation rates. 

Calculated here T ee for J/ip, ip' = ?/>(3686), and 
ip = ijj(377Q) (with the mixing angle 9 — 11°) agree with 
experiment with accuracy < 5% and this allows us to ex- 
tract some important factors from the di-electron widths. 
The essential fact is that correct ratio of the branching 
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0.24, appears 
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to be two times larger then in the "12% rule" mostly 
because different a s describe corresponding annihilation 
rates. 

The unclear situation still remains with two-photon 
width of rj' c , because its value can depend on possi- 
ble influence of virtual decay channel DD* and possi- 
bly other channels [18[ . In closed-channel approximation 
r 71 (r]' c ) = 3 keV is obtained if the same a s = 0.24, which 
provides correct number for T(r)' c — > gg) = 11.2 MeV, 
is taken. This two-photon width is larger than in the 
CLEO experiment [191 ] , where under assumption (uncon- 
firmed ) that B{r\ c -> KKn) = B(r)' c ^> KKw), the value 
r 77 (?7^) = 1.3 ± 0.6 keV has been reported. However, if 
via the DD* channel the mixing of the 3 1 5o and 2 1 Sq 
states occurs, then even with small ^Sq contribution to 
the w.f. of r]' c ( 4% to the norm) its two-photon width is 
becoming essentially smaller, T 77 < 1.9 keV. 



II. VECTOR AND PSEUDOSCALAR DECAY 
CONSTANTS 

The decay constants /y and fp are calculated here 
with the use of the analytic expressions, derived in @. 
To obtain these expressions the functional integral rep- 
resentation for the correlator of the currents (in V and 
P channels) is used and on the final stage this correla- 
tor is expanded in the complete set of the eigenfunctions 
(e.f.) of the relativistic string Hamiltonian (RSH) [20| . 
[2lj . As the first step we use here this RSH to calculate 
charmonium spectrum and define relativistic corrections 
to the decay constants: 
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Here the relativistic factors £p, £,v, refereing to the P and 
V channels, are different and given by the expressions: 
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2w 2 



(3) 

The values of u) and the wave functions (w.f.) at the ori- 
gin are given in Appendix. In fTJ) and i? n (0) (n=l,2) 
refers to the physical radial w.f. at the origin for J/ip, 
ip' = VK3686), and R D (0) is the w.f. of i/>" (3770), i.e. the 
S-D mixing (with the mixing angle 6 — 11°) is taken into 
account. In Appendix for pure 2S and ID states their 
w.f. at the origin are denoted as i?2(0), Rd(0). 

The characteristic feature of RSH is that it contains a 
minimal number of fundamental parameters: the string 
tension cr, Kqcd for n/ = 4, and the pole quark mass m c . 
The string tension is taken from the analysis of the Regge 



trajectories of light-light mesons [2l| and the spectra of 
heavy- light mesons [y] , where the preferable value is a = 
0.180 GeV 2 . The pole mass of c quark is now known 
with rather good accuracy, m c = (1.40 ± 0.05) GeV [22| 
and for nf — A the QCD constant Aj^g = 255(5) MeV is 
taken here. 

The RSH for a meson can be presented as in (23l |: 



UJ 



(4) 



This Hamiltonian H is unperturbed part of general 
Hamiltonian. 



H = Hn + AH, 



(5) 



where AH = VsD{r) + AV str + Vse includes the spin- 
depended part Vsd = Vss + Vls + Vr, the string correc- 
tion V s tr, and the self-energy term Vse, which are consid- 
ered as a perturbation. Notice that in heavy quarkonia 
the string correction (as well as the self-energy term) is 
always small, | V s t r |< 5 MeV, and can be neglected. 
The Hamiltonian Hq has an advantage as compared to 
the spinless Salpeter equation (SSE), since its w.f. at 
the origin for L — is a regular function while the S- 
wave w.f. of SSE diverges at small r [24[ and has to be 
regularized. 

In einbein approximation the spin-averaged mass M u l 
can be presented as: 



M nL = U>r, 



m c 



E nL {uj c ) + Ase, 



(6) 



where the e.v. E n L are the solutions of the so-called 
einbein equation: 



V 1 

— + Vo(r) 



tfinL(r) = E nL ipnL- 



(7) 



Notice that the mass formula ((6|) does not contain 
an arbitrary (fitting) constant Cq. It is essential that 
in einbein approximation for a given state the values, 
M n L and ^> n h , are defined from the extrcmum condition: 
= 0, which provides the accuracy ~ 5 - 7% [23|. 
For the mass ([6]) this extremum condition gives rise to 
the relation: 



OEnL 
d0J n T 



(8) 



Then for a given nL state the dynamical mass u n L and 
M n L are calculated. Notice that u) n L, being the kinetic 
energy of a quark, plays the role of the constituent quark 
mass which slightly differs for the states with different 
quantum numbers. 
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TABLE I: The charmonium spin-averaged masses M n z,(in GeV) with the parameters (m — 1.41 GeV, cxo = 0.18 GeV 2 , 
A SE = -22 MeV, A B (n/ = 4) = 360 MeV, M B = 1.0 GeV)*. 



State 




M nL 


M nL 
a — o(r) 


Experiment 

Mcog 


15 


1.592 


3.068 


3.066 


3.068 


25' 


1.652 


3.678 


3.670 


3.674(1) 


IP 


1.618 


3.488 


3.484 


3.525(1) 






Above DD threshold 




ID 


1.650 


3.787 


3.779 


3.771(2) t 




l.Doo 


3.954 


3.940 




35 


1.712 


4.116 


4.093 


4.040 | 


2D 


1.715 


4.189 


4.165 


4.153(3) T 


3P 


1.742 


4.338 


4.299 




45 


1.772 


4.482 


4.424 


4.421(4) | 


3D 


1.772 


4.537 


4.475 




55 


1.826 


4.806 


4.707 


~4.664(16) 



The symbols (f, J,, ~) mean that not all members of a multiplet are 
measured and therefore the center of gravity cannot be accurately 

defined. 



The masses given in Table 1 are calculated with 
the static potential which contains perturbative gluon- 
exchange (GE) term and nonperturbative confining term: 



t r / x 4a B (r) 

Vo(r) = —z + or, 

S r 



(9) 



where the vector coupling in coordinate space a sir) is 
defined as in 12511: 



otB{r) = - 

7T 



sin(gr) 
dq — a B (q), 



ot B {q) 



o 

47T 
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(3i In i B 



/3o 2 



tB 



(10) 



with t B = In 



q 2 +Mj 



Here M B (a,A B ) = (1.00 ± 0.05) 



GeV is so-called background mass [25| . and As(n/) 
can be expressed through A-jjg'- in 2-loop approxima- 
tion A B {n f = 4) = 0.360(10) MeV corresponds to the 
Ajys = 0.254(7) GeV and in this case the freezing value 
a C rit = 0.547. 

Although here we consider low- lying states, to repre- 
sent gross features of the charmonium spectrum and the 
position of higher levels we take into account flattening 
of linear confining potential. This phenomenon occurs 
due to creation of virtual light-quark pairs plj and it is 
becoming essential for higher levels, in particular, for the 
mass of the S 3 ^! state. The flattening of linear potential 



is defined by the analytic function for which the form and 
parameters are taken just the same as in [2l[ , where they 
have been extracted from the light meson radial Regge 
trajectories. The origin of flattening comes from the vir- 
tual qq pairs creation on the surface inside the Wilson 
loop (W(C)), having large size, and due to these virtual 
loops the string tension (as well as the surface) is be- 
coming smaller and dependent on the QQ separations r. 
This potential provides a correlated mass shift down of 
all radial excitations with n > 3. 

From Table 1 one can also see that in the potential with 
er(r)r the masses M(3S),M(4S), M(5S) are shifted down 
by - 20 MeV, ~ 60 MeV, and ~ 100 MeV, respectively, 
and turn out to be close to the experimental values with 
the exception of the ■0(4040), which is strongly affected 
by the nearby S wave threshold. Since the hyperfine 
splitting of the 5 3 Si level is small, < 6 MeV [Hj], its mass 
practically coincides with the centroid mass, M(5S) = 
4.70 GeV, calculated here, and lies close to the mass of 
the Belle resonance Y(4660) 4]. From Table 1 one can see 
that with exception of the IP state all other masses are 
in good agreement with experiment, even for the states 
above the DD threshold. The reason, why only for the 
IP states the centroid mass has smaller value (in einbein 
approximation) needs an additional analysis. 

The V and P decay constants, calculated with the use 
of (P) and ([!]), are given in Tables 2 and 3 together with 
those from some other papers. 
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TABLE II: The decay constants f v (in MeV) of the J/ip, tp' = -0(3686), i/>" = ^(3770) mesons. 
State This paper BL [9] EFG [11] Wang[12] Experiment 



(Rel) (3 V = 1.0 V = 0.72 
(a s = 0) (as = 0.165) 



J/i> 


483 


545 


551 459(28) 


415(6) 


490(7) 


</>' 


357 


371 


401 364(24) 


302(4) 


356(4) 




115 


318 




96(5) 


113(6) 



TABLE III: The decay constants fp (in MeV) of the r/ c , r/' c , r/ C 2 mesons. 



State This paper 


BL [9] 
(Rel) 


LP [10] 


(3 P = 1.0 
(as = 0) 


Experiment 

a s (rjc) = 0.195 
a.(f£) = 0.25 


Vc 453 


493 


480 


404(57) 


454(64) 


j]' c 336 a) , 267 b) 


260 


303 


189(40) 


213(47) 


Vc2 41 











a > The influence of virtual decay channels is neglected. 
Mixing of 2 1 5'o and S 1 So states (6 — 11°) is taken into account. 



It is worth pointing out that the "experimental" fv 
and fp, extracted from the experimental di-electron and 
two- photon widths, depend on chosen values of radiative 
corrections, 0y and (3 p. Therefore in Tables 1 and 2 
we give two variants of "experimental" decay constants, 
which correspond to /3y = 0.72 (a s = 0.165) and 1.0, 
and (3 P = 0.79 (a s = 0.195) and 1.0 (a s = 0). It is also 
important that 7^(2 So) and ^(l 1 ^) cannot be mixed 
and therefore their w.f. at the origin are defined by the 
w.f. i?2s(0) and i?u(0) for pure 2S and ID states (see 
Appendix). From Table 3 one can see that for /3y = 
0.72 our decay constants fv are in good agreement with 
experiment. Note that for tp" calculated constant fv is 
almost three times smaller than in 191. 



III. DI-ELECTRON WIDTHS 



origin and also several m.e.). In Table 4 calculated widths 
are presented together with some theoretical predictions 
and experimental data. 

Since the coupling in (fT2"|) appears to be the same for 
J/ip, ip' ', tp" , it is of interest to compare their ratios 
(where the radiative corrections are cancelled): R^i = 
r ee (f)/T ee (J/^), iy = r ee (0")/r ee (J/0) with exper- 
imental numbers which turn out to be very close to each 
other: 



(ity)t/i = 0.46, (iV W = 0.45 ± 0.02, 
(iV'k = 0-046, (JV W = 0.044 ± 0.03. 



(13) 



Having such agreement with experiment we expect 
that our w.f. at the origin are defined with good ac- 
curacy, in particular, the following numbers are obtained 
for the factors : 



The leptonic width of a vector state in heavy quarkonia 
is expressed via the decay constant [f|,[26j]: 



i 



3My 

16 

3^- 



(11) 
(12) 



Best description of di-electron widths is obtained here 
taking in (fTT|) a s = 0.165, or [iv = 0.72, for which 
r ee (J/0) = 5.41 keV, r ee (V/) = 2.47 keV, andr ee (^") = 
0.248 keV are obtained (see Table 4). In our calculations 
of the di-electron widths of ?/>'(3686) and V"(3770) the 
S-D mixing with the mixing angle 9 = 11° is taken into 
account (In Appendix we give their physical w.f. at the 



i(0) 



A/1 



= 0.085Gey, 



K'(Q)| 2 - 



= 0.0040Gey, 



Km 



= 0.040GeV, (14) 



= 0.0027Gey, (15) 



We estimate the accuracy for these factors as < 10% and 
later use them to define three- and two-gluon annihilation 
rates for J/ip, tp', i] c , r( c . 

Notice that relativistic corrections decrease the decay 
constants fv and fp and provide better agreement with 
experiment. However, for hadronic decays this type of 
relativistic corrections is not still calculated in FCM, and 
therefore the accuracy of calculated hadronic widths is 
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TABLE IV: The di-electron widths (in keV) (with a s = 0.165 
and the mixing angle 9 = 11°). 

State This paper BL[27] EFG [11] Experiment 



J/i/> 



5.41 
2.47 
0.248 



12.13 
5.03 
0.056 



5.4 
2.4 



5.55 ±0.16 
2.48 ± 0.06 
0.242 



+0.027 
0.024 



( ^77 (Xcq ) \ 



0.41 ±0.02. 



For Xco the width is by 10 % larger than the experimental 
one, while for Xc2 the absolute value and the ratios: 



worse than for the decay constants fy and fp. Con- 
cluding this Section, we would like to underline that our 
calculations give the ratio of the branching fractions for 



/ r (x ) \ _ 07g (x ) \ 



ex p 



0.076 ±0.002 

coincide with the experimental values. Different two- 
photon widths are given in Table 5. 



Re 



geeC^Q _ r t ot(J/lp) \R^(0)\ 2 M 2 (J/j>)Z v {2S) 

= (12.6 ±0.03)%, 



(16) 

which is in good agreement with experimental number, 
or the so-called "12% rule", R ee (exp) = (12.3 ± 0.03)%. 



IV. TWO-PHOTON WIDTHS OF 

r7 c ,77 c , Xc (l 3 Po),Xc(l 3 P 2 ) 

The two-photon widths of rj c (lS), rf c (2S) can be ex- 
pressed via the decay constants fp: 



TABLE V: The two-photon widths (in keV) with a s = 0.195 
for 77c, Xco, Xc2 and with a s = 0.24 for r)' c , r) C 2- 



State This paper BL [9] EFG [11] KLW [12] Experiment 



r? c 7.22 
r)' c 3.0 a) , 1.9 b) 
r? c2 0.042 
Xc0 (l 3 Po) 3.31 
Xc 2 (l 3 P 2 ) 0.54 



4.18 
2.59 
1.21 
3.28 



5.5 7.14(95) 7.1±2.7[22J 

1.8 4.44(48) 1.3 ±0.6 [18] 

2.9 3.78 2.90 ± 0.43 
0.52 0.53 ± 0.06 



T 77 (P) = -Mffr^ 

Pp = 1 5 a s . 

3tt 



(17) 
(18) 



The two-photon widths of the scalar \ c q and tensor \c2 
mesons can be derived in FCM the same manner as it 
has been done for the V, P decay constants in [6j with 
the following result: 



f 2 



S(T) 



3\R' P (0)f 



S(T)U P 



irM 



(19) 



Here in (I19p instead of the quark mass uiq (in nonrela- 
tivistic limit) the kinetic energy ujp enters. Then with the 
QCD corrections the two-photon widths of the P— wave 
mesons are defined as [7j: 



-T77(Xco) 
r 77(Xc 2 ) 



108a 2 e4|j^(Q)| ; 
Mjujp 

144a 2 ef|^(0)r 



28 - 3tt 2 



Oh 
7T 



1 - 



16 as 

3 7T 



Q) 

(21) 



In our analysis we use a s — 0.195 for r/ c , XcOi an d 
Xc2 to obtain good numbers for the following ratios of 
two-photon widths (in which relativistic factors £,p(r] c ) = 
0.785, £p(r)' c ) — Cp(^c2) = 0.73 are used, see Appendix): 



a) See the footnote a) to Table 3. 
6) See the footnote b) to Table 3. 



We would like to stress here that calculated two-photon 
width of rj' c (with a s = 0.24) is larger than in the CLEO 
experiment [19J , nevertheless with the same a s we have 
obtained hadronic width T(rj' c — » qg) = 11.0 MeV in 
agreement with T tot = 14 ± 7 MeV |22T ]. 

One cannot exclude that the w.f. at the origin of rj' c is 
affected by the virtual DD* decay channel which lies only 
by 130 MeV higher than rf c . Then via this channel the 
w.f. of rj' c can be mixed with ^Sq state (which is now 
often identified with A" (3940) For example, the 

20% admixture in the w.f. of rj' c (or the 4% contribution 
to the norm) gives rise to the two-photon r 77 (?y^) = 1.9 
keV. 

Recently two- pho ton width of r\ c has been calculated 
in lattice QCD |30l], [3lT ]. Such calculations from first 
pri ncip les are very important for the theory, however, 
in df rather small number, r 77 (r7 c ) = 2.65(26)(80)(53) 
keV and corresponding fp{rj c ) — 373 MeV are obtained 
in quenched approximation. Meanwhile in [31( the calcu- 
lations in lattice QCD with exact chiral symmetry, where 
heavy quarks are treated as the Dirac fcrmions, the larger 
decay constant, fp(r) c ) — 438(11) MeV, is obtained and 
this number is in agreement with our result, fp(T] c ) = 453 
MeV. 
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THREE- AND TWO-GLUON 
ANNIHILATION RATES 



We would like to notice also that in our calculations 
the ratio, 



The QCD corrections to di-electron widths have ap- 
peared to be < 30%, but they are even more important 
for some hadronic decays. From Q, @|, [3 we know 
the widths of the three-gluon annihilation and the decay 
into rygg for the vector mesons, as well as for two-gluon 
annihilation of P, S, T mesons: 



T(V- 
T(V - 
T{P 

r (xc 
r (xc 2 



999) 
199) 

-* 99) 
99) 

■* 99) 



40(tt 2 - 9)a 3 s 



\Rn(0)\ 2 



32(n 2 



9)a2sa \Rn(0)\ 2 



1-3.7— (24) 

7T ' 



6.7^5) 

7T 



3M P 
24a 2 . 

32a? 



hM T u? p 



|i?„(0)| 2 (l + 4.8^ 



|^(0)| 2 1 



9.5^ 

7T 



\R 



p(o)| 2 ( 



1 - 2.2- 



(26) 
(27) 
(28) 



It can be easily shown that with a s = 0.165, as for 
considered above di-electron widths, the three-gluon an- 
nihilation width of J/ip is smaller than in experiment 
being equal 42 keV. The best fit to this annihilation rate 
is obtained taking a s = 0.187 (practiclly the same as in 
two-photon width), for which 

T(J/iP^ggg) = 59.5 keV, 
T(J/ip -> 19 g) = 5.7 keV. 

Then the sum of these annihilation widths , being equal 
the width T(J/tp — » light hadrons) = 65.2 keV, is in 
good agreement with the experimental number B( J/ip — > 
light hadrons) = (69 ± 3)% 0, [13 or T(J/ip -> 
light hadrons) exp — 64 ± 3 keV. 

On the other hand this value, a s — 0.187, is not suf- 
ficient to provide correct number for the three-gluon an- 
nihilation rate of ■0'(3686) and in this case the best fit is 
obtained taking a s = 0.238: 

T^'^ggg) = 52.7 keV, 
r(V' -> 735) = 3.5 keV, 

with their sum T(ip' — > light hadrons) = 56.2 keV, in 
good agreement with the experimental branching fraction 
B(ip' -> light hadrons) = (16.9±3)% 0, 0, or T(ip' -> 
light hadrons) ea; p = 57±10 keV. This fact means that the 
ratio 



Vto 



d _ g(7/i'->light hadrons) _ 
LH ~ B(j/v-^light hadrons) — 
{J/ip) r(v»'->light hadrons) _ q 24 
r(j/V'^light hadrons) — 



(29) 



is in agreement with experiment and two times larger 
than the ratio R ee (1 mostly because different a s are 
taken for J/ip and ip' . 



= fW -> light hadrons) = ^ Q ± Q ^ 

r ee (j/ip) 

agrees with experimental number, i? eKp = 11.6 ± 0.03, 
and this fact justifies our choice of a s = 0.187 in radiative 
correction for J/ip. 



VI. CONCLUSIONS 

We have shown that di-electron widths of J/ip, ip 1 , ip 
describe experimental data with the accuracy better 5%, 
if the same a s — 0.165 is taken in the QCD radiative 
corrections. This fact can be considered as a test of the 
method used here. 

For i] c (lS), Xco, Xc2 the larger effective coupling, 
a s = 0.195, is needed to fit experimental numbers for 
two-photon widths. Also close value of a s = 0.187 pro- 
vides good description of the annihilation widths for 
the J/ip — ► ggg and J/ip — > jgg processes, so that 
B(J/ip — > light hadrons) = 70% is obtained (its value 
in experiment is (69 ± 3)% [3, [l7|). 

However, to describe experimental branching fraction, 
B{ip' -> light hadrons) = 16.7 ± 3.0% [II], 0, larger 
a s = 0.238 is needed, for which we obtain T(ip' — > <?g<?) = 
52.7 keV and r(-0' — * 793) = 3.5 keV, so that their sum 
is equal T(?// — -> ZigM hadrons) — 56.2 keV, which cor- 
responds to experimental branching fraction. Then the 
ratio of the branching fractions: 



-Riff — 



— > light hadrons) 
B(J/ip — > light hadrons) 



= 0.24 ±0.01, 



appears to be in good agreement with experimental num- 
ber, Rlh(&xp) = 0.24±0.04, being two times larger than 
the ratio R ee (fT6|) . 

Nevertheless for the same a s — 0.24 calculated here 
two-photon width of rf c , T- y y(r) / c ) — 3.0 keV, turns out to 
be larger than in the CLEO experiment [lJ|. From our 
point of view the problem of this width can be solved tak- 
ing into account the influence of the DD* decay channel. 

Thus our analysis shows that there is no an universal 
effective strong coupling which allows to describe differ- 
ent annihilation processes for all low-lying charmonium 
states. Such universal description is possible only for di- 
electron widths. 
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APPENDIX: THE MATRIX ELEMENTS AND 
WAVE FUNCTION AT THE ORIGIN OF 
LOW-LYING CHARMONIUM STATES 

Here we use RSH to calculate the m.e., like u>(nL), < 
p 1 >, and also the w. f. at the origin. For pure 2S state 
we denote its w.f. as R 2S (0) = 0.767 GeV 3 / 2 ; for the 
P-wave states R P (Q) = R' p (0)/ujp = 0.183 GeV 3 / 2 with 
R' P (0) = 0.297 GeV 5 / 2 , and for the ID state R D (0) = 
|^2$£ = 0.0942 GeV 3/2 with i?^(0) = 0.145 GeV 7/2 . 

The values of Lo(nL) are given in Table 6. 

The physical w.f. of ip' and ip", calculated here, take 
into account the mixing angle = 11°, and these w.f. at 
the origin are given in Table 6 (all needed parameters are 
given in Section 2). 

For ip 1 and ip the S-D mixing can occur due to tensor 
forces and coupling to open DD channel. Then the w. f. 
at the origin of ip' and ip are defined as in [28| : 

i^,(0) =cos0J2 2 ,s(O) =—sin0i?£(O), (30) 



iV'(0) = sin 0R*s(0) + i - cosPi%(0). (31) 
2v2ur 

Our calculations give Pp(0) = |fl ^ ( ° )|2 = 0.183 
GeV 3 / 2 ; the V" w.f. at the origin, i? D (0) = 0.238 
GeV 3 / 2 (w(l-D) = 1.65 GeV), and for the ip' meson 
iV(0) = 0.735 (see Table 6). 

TABLE VI: The dynamical masses ui n L , the radial wave func- 
tions at the origin for J/ip ,ip' , ip , and \ c j (in GeV). 



State uj nL , GeV R nL (0) <f>jv jp 



IS 


1.59 


0.905 


0.541 


0.929 0.785 


2S 


1.65 


0.735 


0.722 


0.910 0.733 


IP 


1.62 


0.183 


0.619 


1.0 1.0 


ID 


1.65 


0.238 


0.721 


0.911 0.733 
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